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Abstract 

A forward dispersion calculation is implemented for the spin polarizabilities 
71, • • • , 74 of the proton and the neutron. These polarizabilities are related to 
the spin structure of the nucleon at low energies and are structure-constants 
of the Compton scattering amplitude at 0(w 3 ). In the absence of a direct 
experimental measurement of these quantities, a dispersion calculation serves 
the purpose of constraining the model building, and of comparing with recent 
calculations in heavy baryon chiral perturbation theory. 
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1. Introduction. 

The response of the internal degrees of freedom of the nucleon to an external electro- 
magnetic field can be parametrized in terms of the structure dependent polarizabilities. The 
classic process for studying such quantities is Compton scattering at energies below the reso- 
nance region. In an expansion of the Compton scattering amplitude in powers of the incident 
photon energy u, the 0(1) and 0(u) terms depend only on the mass M, the electric charge 
e and the anomalous magnetic moment k of the nucleon. Therefore, no information on the 
excitation of the internal degrees of freedom can be obtained up to O(ou); it is only at O{oj 2 ) 
that the amplitude becomes sensitive to the internal excitation of the nucleon. At 0(u 2 ), the 
amplitude is parametrized in terms of the electric (a) and magnetic {(3) polarizabilities [I], 
which describe the deformation of the system in the presence of a static electric and magnetic 
field. At 0(u 3 ), as demonstrated recently by Ragusa 0, four new polarizabilities appear, 
7ij 72)73; an d 74, the "spin polarizabilities". As explained by Ragusa, these can be inter- 
preted as the response to the external fields of a magnet that has an internal structure. For 
related previous work on the 0(u> 3 ) spin structure of the Compton scattering, see Ref. |§. 

Contrary to a and /3, none of the ji's has been measured experimentally, although there 
has been a proposal for that purpose . Only the spin-dependent polarizability for forward 
scattering involving three of the spin polarizabilities, 7 = 7i — 72 — 274, has been experi- 
mentally constrained through an analysis involving photoproduction multipoles || , with the 
result 7 = —1.3 (here and thereafter we quote all polarizabilities in units of 10 _4 fm 4 ). An 
earlier multipole analysis || obtained 7 = — 1.0 . 

On the theoretical side, there are recent calculations of these quantities using chiral per- 
turbation theory (ChPT). In Ref. f7j, 7 was evaluated in a 0(p 3 ) calculation in heavy baryon 
ChPT (HBChPT), with the result: 



' 4vr 247r 2 F2m2 ' ' v ' 

where F n = 92.4 MeV, m n = 137 MeV, and g^ = 1.26. This prediction is in clear disagree- 
ment with the multipole analysis. Sub-leading effects evaluated within a relativistic ChPT 
calculation at the one-loop level diminish the discrepancy, yielding <-y 1 ~' oop = 2.2. More- 
over, inclusion of the A(1232) resonance through higher order contact terms within the same 
approach [7| revealed large effects on 7, with opposite sign to the one-loop result, j A = —3.7. 
The net result in the relativistic calculation, 

7 = +2.2 (1-loop) -3.7(A) = -1.5, (2) 
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is very close to the one found with the multipole analysis mentioned above. 

The individual spin polarizabilities were calculated in Ref. || within the same HBChPT 
approach as in Ref. 0j. An interesting result of this calculation is that the values of 71,73 
and 74 are completely dominated by the contribution of the t-channel 7r°-exchange (Wess- 
Zumino-Witten term), while this contribution cancels in the expression for 7. Specifically, 
the expressions obtained for the 7j's are 

_ e 2 g A f 9a\ _ e 2 g A / g A 

71 47747^774 V 6 J ' 72 47r47r 2 F2 m 2 V 12 

e 2 gA /l . 9a\ 

73 = 9P9 9 n + 7T7 » 74 = -73- (3) 



4tt 47r 2 F2m2 V 2 24 

Very recently, in Ref. || the A(1232) was introduced as an explicit degree of freedom 
in HBChPT. In this approach a new dimensionful parameter enters, A = — M, and a 
consistent HBChPT formalism can be set up according to an 0(e n ) power counting scheme, 
where e denotes a small scale (a soft momentum, or A). In this formalism the previous 
0(p n ) HBChPT results || are exactly reproduced, and the terms that come from treating the 
A(1232) as an explicit degree of freedom can be clearly identified. The value for 7 obtained 
to C(e 3 ) is 



e 2 



- f 9a 4& i 1 Ana 

7 4vr I 24^^27772 9M 2 A 2 54tt 2 F2 
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A 2 + 2m2 3m 2 A • 

In it 



(4) 



_(A 2 -m2)2 (A 2 -7772)5/2 

where R = A/m n + [(A/m^) 2 — 1] 1//2 , ^ata is the pion-nucleon-A(1232) coupling constant, 
and bi is a finite 0(e 2 ) counterterm. The terms in Eq. (^) correspond, in order, to nucleon- 
loop (N-loop), delta-pole (A-pole) and delta-loop (A-loop). As remarked in Ref. |J, both 
gnNA and bi are not well constrained by the experiment. However, in using the values 
g^NA = 1.5 ± 0.2 and 61 = -(2.5 ± 0.35) in Eq. (§, Hemmert et al. § obtained 



7 = 4.6 (N-loop) - 4.0 (A - pole) - 0.4 (A - loop) = +0.2. (5) 

As in the calculation of Ref. 0, the contribution of the A (1232) is large and negative, but 
cannot account for the value 7 = —1.3 found from the multipole analysis. Unfortunately, 
also the values for the electric and magnetic polarizabilities, a and j3, are larger than the 
experimental data. 

More recently, Holstein et al. [|ll|] |L2| have argued that the values of g n NA and b\ to be 



used should be determined, for consistency, within the same "small scale expansion" of the 
HBChPT calculation, and not within a relativistic Born model. In such a case, the values for 
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these quantities are approximately 50% smaller than the ones used above. The consequence 
of this is that a and (3 decrease significantly, but on the other hand, the value for 7 increases 
significantly: 

7 = 4.6 (N-loop) - 2.4 (A - pole) - 0.2 (A - loop) = +2.0, (6) 

in contradiction with the multipole analysis of Ref . H . The effect of the A pole is also quite 
large for some of the 7,'s. In Table I we show the results for the 7;'s using for g n NA and h\ 
the values suggested by Holstein et al. |TT| |T2"]. 



As mentioned above, there are no direct measurements of the spin polarizabilities. Nev- 
ertheless, as in the case of 7 || ||, an estimate for the 7j's can be obtained using the experi- 
mentally determined amplitudes for pion photoproduction in a dispersion integral. However, 
there are additional contributions from t-channel processes that are not well constrained 
experimentally, corresponding to possibly large high-energy contributions to the dispersion 
integrals, similar to the effects of the Wess-Zumino-Witten term in the case of ChPT. In this 
sense, only the low energy part of the Compton process, i.e. contributions of the pion cloud 
and the low-lying resonances, can be at present constrained by a dispersion calculation. 

In this Letter we present the results of a dispersion calculation for the spin polarizabilities. 
In order to minimize uncertainties with respect to the high-energy behavior of the Compton 



amplitude, we follow the approach of L'vov et al. |I3|. In this approach one uses a finite- 
size contour in the complex plane rather than improving the convergence of the integral by 
subtractions. The contributions arising from the contour are expressed in terms of low-mass 
meson exchanges in the t-channel, as will be discussed in the next section. We use pion 
photoproduction multipoles from a recent analysis based on fixed t dispersion relations and 
unitarity by Hanstein, Drechsel and Tiator (HDT) [|l4j . This analysis is limited to a maximum 
photon laboratory energy of the order of 500 MeV. In order to investigate the sensitivity of our 
results to the photoproduction amplitudes at higher energies, we also employ the multipole 
analysis VPI-SP97K of the Scattering Analysis Interactive Dial-in program, SAID ||15|| . In 
Section 3 we present our numerical results and compare them with the results of ChPT. Some 
conclusions and perspectives are presented in Section 4. 
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2. Dispersion relations for the spin polarizabilities. 

In the cm. system the T-matrix for Compton scattering on the nucleoli can be written, 
in the Coulomb gauge, in terms of six amplitudes that are functions of the energy of the 
incident photon, u, and the scattering angle, 9, as [[RJ: 

T = xt> {e'*.eA 1 (o;,0) + e'*-ke-k' A 2 (u,0) + i<x-(e'*xe) A 3 (oo, 9) + ia-{k' Xk)e'*-e A a (lu, 9) 
+ i(T'[(e'*Xk)e-k' - (exk')e'*-k] A 5 (u,9) + ia-[{^*xk')e-k' - {exk)e'*-k}A 6 {u,9)} X s, (7) 

where k, e (k', e') are the direction and the polarization vector of the incident (final) photon, 
Xs {Xs') is the initial (final) nucleon spinor, and cr is the Pauli spin matrix. For convenience, 
we write the amplitudes Aj as a sum of the Born and non-Born contributions, A^ = Af orn + 
Af B . Ragusa's spin polarizabilities can be expressed in terms of the non-Born amplitudes 
^ B 's by § 



i i d 3 



An 6 dui 3 

i i a 3 



Af(u, 0) + A^(u, 0) + 2Af(u, 0) 



nB/ 



1 1 d 3 



73 = ^—7 



47T 6 did 3 



uj=0 

1 1 d 3 



72 = ^7 



47r 6 duj 3 



uj=0 



uj=0 



74 = ^7 



47T 6 did 3 



Af(u,0) 



uj=0 



On the other hand, the scattering amplitude can be written in terms of the Hearn-Leader 



(HL) |T7| Lorentz-scalar amplitudes Tj(z/, t), 



1 • • • 6, as 



i=l 



u s {p) e L 



(9) 



where e M and e v are the polarization four- vectors of the nucleon, u s (p) and u s (p) are Dirac 
spinors, and P are tensors that depend on Dirac 7 matrices and the initial and final momenta 
of the photon and the proton. From arguments based on the asymptotic behavior of Regge 
trajectories, some of the HL amplitudes appear to have a bad convergence for fixed t and large 
v. The traditional approach in dispersion theory is to implement a subtraction at threshold 
for these amplitudes with a bad asymptotic behavior, at the cost of introducing subtraction 
functions that in most cases are unknown. However, there is an additional problem with the 
HL amplitudes Ti(v,t), in that they have to fulfil kinematic constraints because the tensors 



P develop singularities at forward and backward scattering angles |L8| . It is possible to avoid 



kinematic constraints and singularities by using appropriate combinations of the amplitudes, 



as for example the ones of Bardeen and Tung ||18|| . One natural way to proceed is the one 
followed by Pfeil, Rollnik and Stankowski |19|| , who used a partial wave decomposition of the 
HL amplitudes. While the constraints can be automatically satisfied within this procedure, 
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it has the disadvantage of introducing a possible violation of the s — u crossing symmetry. 
Another approach has been introduced by L'vov et. al |13|], who expresses the dispersion 
contribution as the sum of a finite-range integral from threshold z/thr up to a maximum value 
^max, and the contributions from higher energies in terms of t-channel poles. Thus, the real 
parts of the amplitudes are written as a sum of the real pole Born terms, the contribution 
of the finite-range dispersion integral, and an asymptotic t-channel contribution (we use the 
conventions and notations of Ref. ITO 



ReAi(u,t) = A\ 



Born, 



V, t) + -P 

7T 



v' dv' 



yl2 _ 



(10) 



The amplitudes A^u, t) are appropriate combinations of the HL amplitudes, and are free 
of kinematic constraints and singularities. Two of these six amplitudes, Ai and A 2 , appear 
to have a bad convergence behavior for high v, at fixed t. The t-channel contribution to 
A*f(t) and Af(t) are modeled by o and 7r° exchanges, respectively, and for the remaining 
amplitudes the t-channel contributions seem to give negligible contribution at low energies. 

In order to evaluate the spin polarizabilities in terms of the dispersion integral, one needs 
the relation of the amplitudes Ai(u,8) of Eq. ([?]) to the Lorentz-scalar amplitudes Ai(v, t). 
The relation is easily obtained comparing Eq. (^) with the low energy expansion of the 
Compton scattering amplitude written in terms of the non-Born contributions to the Lorentz- 
scalars Ai(v, t), with the result |2(| [2]J 





1 1 r 


7i = 


4ttM - 




1 1 r 


72 = 


4^M - 




1 1 r 


73 = 


+ 8ttM - 




1 1 r 


74 = 


~8%M - 



Af (0,0) -Af (0,0) + Af (0,0) 
Af (0,0)- Af (0,0) 
Af (0,0) -Af (0,0) -Af (0,0) 
Af(0, 0) + Af(0, 0) + 2Af (0, 0) - Af(0, 0) 



From these expressions one sees that the polarizabilities 71, 73 and 74 depend on the amplitude 
A 2 . As we discussed previously, this amplitude has a bad asymptotic behavior for fixed t and 
large v and not well-constrained t-channel contributions. One way to avoid such uncertainties 
is to consider appropriate combinations of the 7i's to which A 2 does not contribute. In 
addition to the combination that leads to 7, there are more of such combinations, which will 
be discussed in the next section. 
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3. Numerical results. 

In the dispersion integral Eq. (|H]) we saturate the imaginary parts by one-pion photo- 
production amplitudes. We neglect intermediate states with more pions and heavier mesons. 
The contributions of such states to 7 have been estimated in Ref. || to be relatively small. 
As stated above, we employ two sets of photoproduction amplitudes, the ones obtained with 
a fixed t dispersion relations (HDT) and the ones from the SAID program. Our results using 
the dispersion relations (DR) for the spin polarizabilities of the proton and the neutron are 
presented in Tables I, II, and III. 

The results in Table I are obtained by evaluating the integral in Eq. fliTf ) with the 
HDT @ multipoles up to a z/ max = E™ ax + t/4M, with t = and E™ 8 * = 500 MeV. 



Tables II and III also show the results for the integration with the VPI-SP97K SAID multi- 
poles up to _E™ ax = 500 and 1500 MeV, which are denoted by SAIDi and SAID2, respectively. 
For the t-channel contribution A^(t) we use the parametrization given in Ref. [13| . 



In Table I we show the separate contributions from tt° exchange and the dispersion inte- 
grals ( "excitation" ) to the polarizabilities. The contribution of the 7r° exchange to Eq. (|T0"D is 
practically identical to the Wess-Zumino-Witten term of the HBChPT calculation. Because 
of the huge contribution from that term, our global result ("sum") is quite similar to the 



prediction of HBChPT ||12|| . However, the contributions beyond the anomaly show significant 
differences. In Table II we present the results for these excitation contributions to the spin 
polarizabilities using the HDT and SAID multipoles. One notices that the HDT and SAIDx 
results are fairly similar for all the 7j's. On the other hand, extending the upper limit of the 
integration of the SAID multipoles from 500 MeV to 1500 MeV changes the results for Ji\ 
^2 P \ an d 74 by a factor of the order of 25%, and changes the sign and magnitude of 73 . 
The values of the 7j- n ^'s are more stable with respect to the change of the upper limit of the 
itegration. We come back to this point shortly ahead. As already seen in Table I, the results 
from DR are at variance with the predictions of HBChPT in most cases. 

As we discussed above, there exist combinations of the 7j's that do not depend on the 
amplitude A 2 and thus are not affected by the badly known high-energy contributions. One 
of these combinations is 7 = 71 — 72 — 274, and three other combinations are: 

713 = 71 +273, 714 = 71 - 274, 734 = 73 + 74- (12) 

Note that only two of these combinations are independent, since 713 — 714 = 2734. The results 
for 7, 713, and 714 are presented in Table III. 
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Our results for and 7^ using the SAID multipoles up to 1500 MeV are practically 
identical to the values obtained for these quantities in Ref. || using the VPI-FA93 analysis 
(see Ref. § for details). We note that the values of 7 (p) and 7 (n) using the VPI-SP97K SAID 
multipoles up to 500 MeV and 1500 MeV differ by less than 10%, i.e. there are no large 
contributions to the dispersion integrals from energies between 500 MeV and 1500 MeV. 
This fact can be understood as due to the damping factor 1/cu 3 in the integrand of the 
expression for 7, 

[<*> a 1/2 - (T3/2 , M „v 

where a%/2 and a 3 / 2 are the total photoabsorption cross sections measured with the photon 
and nucleon polarizations parallel and antiparallel, respectively. On the other hand, we find 
that the value of 7 calculated with the HDT multipoles is a factor of two smaller than the 
value calculated with the SAID multipoles. This can be traced to differences between the 
E 0+ amplitudes in both analyses. In Figs. 1 (a) and (b) we plot the contributions of the 
E 0+ and M 1+ multipoles (which are by far the largest contributors) to the integrand in 
Eq. (0), and in Fig. 1(c) we plot the total integrand. From these figures it becomes obvious 
that the difference comes from the behavior of the E 0+ multipoles close to threshold. As 



remarked in Ref. [22], the SAID multipoles are not meant to be used in the 7r + n threshold 
region. For example, the amplitude E Q+ (n7T + ) at threshold is 24.9 x 10 _3 /m 7r + for SAID and 
28.3 x 10~ 3 /m 7r + for HDT, the latter being much closer to the threshold value predicted by 
ChPT @, 28.4 x 10" 3 /m^+. The differences between the E 0+ amplitudes of the HDT and 
SAID multipoles is of the order of 40% at the peak value of the integrand of 7. As seen in 
Fig. 1, this large difference between the E Q+ amplitudes, and a small difference between the 
M 1+ amplitudes (the SAID value is a little larger than the HDT value), give a net effect of 
50% in the final value for 7. 

An interesting comparison can be made between the DR and HBChPT results by con- 
sidering in separate the leading contributions to 7, which come from the E 0+ , E i+ , and M 1+ 
multipoles. Using the HDT multipoles, we obtain: 

7 = +2.5 (E 0+ ) - 3.0 (E 1+ , M 1+ ) - 0.1 (rest) = -0.6, (14) 

where the last term includes all other partial waves. This shows the strong cancelation 
between s-wave pion loop and A resonance contributions, which is in close correspondence 
with a similarly strong cancelation in the ChPT calculations, Eqs. (|2|), ([]), and (](]). It is 



amusing to note that the first two numbers in Eq. (jT4|) are not very different from the ones 
in Eq. (|2|) obtained in the relativistic ChPT calculation of Ref. |7[], namely ■y 1 ~ loop = 2.2 and 
7 A = -3.7. 

To conclude our discussion on the numerical results, we note that the results in Table III 
clearly show that when the troublesome amplitude A 2 (V, t) is eliminated by taking appropriate 
combinations of the 7i's, the results become very stable with respect to the upper limit of the 
integral. In fact, similar to the case with 7, the integrals for the other combinations of 7j's 
are almost completely saturated with an upper limit of 500 MeV. The large changes in the 
SAIDi and SAID2 predictions observed in Table II can therefore be traced to the presence 
of the function A 2 (i > , t) in the expressions of the 7j's. In fact, inspection of the integrands 
for the Tj's reveals that these still receive sizable contributions from high energies due to 
the function A^iy, t). 

4. Conclusions and Perspectives. 

In this paper we use dispersion relations to estimate the values of the spin polarizabilities 
7i of the nucleon and compare our results with HBChPT calculations. By taking appropriate 
combinations of the 7i's, we are able to minimize uncertainties related to t-channel contri- 
butions. The present dispersion calculation predicts, in accord with previous calculations H 
M, an opposite sign for 7 = 71 — 72 — 274 as compared to the HBChPT calculation. The 



same discrepancy is seen for 714 = 71 — 274. Note that in the HBChPT calculation [12 

71 receives sizable contributions from tt — N loops but no A-pole contributions, whereas in 

72 and 74 there are cancelations between contributions from it — N loops and the A-pole 
graphs ||12|| . Therefore the required change of sign in both 7 and 714 seems nontrivial and 
highly constrained. In this sense, it would be extremely interesting to see the outcome of a 
systematic HBChPT calculation of these quantities to 0(e 4 ). 

Another interesting conclusion of our calculation is that the spin polarizabilities are very 
sensitive to the behavior of the photoproduction multipole E 0+ near threshold. There is 
a large cancelation between the contributions from the Eq + and M\ + multipoles, and a 
very precise determination of both of these is necessary to constrain the values of the spin 
polarizabilities. 

We conclude by remarking that similar to the electric and magnetic polarizabilities a 
and /3, the spin polarizabilities 71 ■ • • 74 contain vital information on the low energy structure 
of the nucleon. The experimental determination of these new polarizabilities is therefore of 
great interest. 
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TABLES 



TABLE I. Separate contributions to the spin polarizabilities from the HBChPT calculation of 
Ref. [12] and the result from the dispersion calculation using the HDT [14] multipoles (all results 
are in units of 10 fm ). 





HBChPT 


DR (HDT) 


li 


wzw 


excitation 


sum 


vT°-exchange 


excitation 


sum 


ip) 

Ti 


-22.0 


+4.4 


-17.6 


-22.5 


+5.1 


-17.4 


ip) 

7 2 





-0.3 


-0.3 





-1.1 


-1.1 




+11.0 


+1.1 


+12.1 


+11.2 


-0.6 


+10.6 




-11.0 


+1.3 


-9.7 


-11.2 


+3.4 


-7.9 


(n) 

7i 


+22.0 


+4.4 


+26.4 


+22.5 


+6.1 


+28.6 


W 
7 2 





-0.3 


-0.3 





-0.8 


-0.8 


in) 

7 3 


-11.0 


+1.1 


-9.9 


-11.2 


-0.6 


-11.8 


in) 
74 


+11.0 


+1.3 


+12.3 


+11.2 


+3.4 


+14.6 



TABLE II. Excitation contribution to the spin polarizabilities from the HBChPT calculation 
of Ref. [10] and the result from the dispersion calculation using the HDT, SAIDi, and SAID2 
multipoles (all results are in units of 10~ 4 fm 4 ). 



7i — excit. 


HBChPT 


DR (HDT) 


DR (SAIDi) 


DR (SAID2) 




+4.4 


+5.1 


+4.3 


+3.5 


7 « 


-0.3 


-1.1 


-1.2 


-1.0 


>) 

73 


+1.1 


-0.6 


-0.5 


+0.1 


ip) 

74 


+1.3 


+3.4 


+3.4 


+2.9 


in) 
7l 


+4.4 


+6.1 


+5.9 


+6.1 


in) 
72 


-0.3 


-0.8 


-1.0 


-0.9 


in) 
73 


+1.1 


-0.6 


-0.6 


-0.6 


in) 
% 


+1.3 


+3.4 


+3.5 


+3.6 
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TABLE III. Combinations of spin polarizabilities that do not depend on the amplitude i) 
(all results are in units of 10 -4 fm 4 ). 



7's 


HBChPT 


DR (HDT) 


SAID! 


SAID 2 


7 (p) 


+2.0 


-0.6 


-1.2 


-1.3 


>) 

7l3 


+6.6 


+3.8 


+3.3 


+3.7 


(p) 

7l4 


+1.8 


-1.7 


-2.4 


-2.3 


7 w 


+2.0 


+0.0 


-0.2 


-0.3 


(n) 
713 


+6.6 


+4.9 


+4.7 


+4.9 


W 
7l4 


+1.8 


-0.7 


-1.2 


-1.1 
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FIG. 1. Contribution of the multipoles Eq + (a) and M\ + (b) to the integrand for 7^ (c) from 
the HDT (solid line) and the SAID (dashed line) analyses. 
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